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Abstract 

A nonuniform extension of the Glauber model on a one-dimensional lat- 
tice with boundaries is investigated. Based on detailed balance, reaction 
rates are proposed for the system. The static behavior of the system is 
investigated. It is shown that there are cases where the system exhibits a 
static phase transition, which is a change of behavior of the static profile 
of the expectation values of the spins near end points. 
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1 Introduction 



The Glauber dynamics was originally proposed to study the relaxation of the 
Ising model near equilibrium states. It is a simple non-equilibrium model of 
interacting spins with spin- flip dynamics. It is also known that there is a re- 
lation between the kinetic Ising model at zero temperature and the diffusion 
annihilation model in one dimension. There is an equivalence between domain 
walls in the Ising model and particles in the diffusion annihilation model. Ki- 
netic generalizations of the Ising model, for example the Glauber model or the 
Kawasaki model, are phenomenological models and have been extensively stud- 
ied [1~6]. Combination of the Glauber and the Kawasaki dynamics has been 
also considered [7-9]. Most studies are focused on completely uniform lattices 
with site-independent rates. Among the simplest generalizations beyond a com- 
pletely uniform system is a lattice with alternating rates. In [10-12], the steady 
state configurational probabilities of an Ising spin chain driven out of equilib- 
rium by a coupling to two heat baths has been investigated. An example is a 
one-dimensional Ising model on a ring, in which the evolution is according to a 
generalization of Glauber rates, such that spins at even (odd) lattice sites ex- 
perience a temperature T e (T D ). In this model the detailed balance is violated. 
The response function to an infinitesimal magnetic field for the Ising-Glauber 
model with arbitrary exchange couplings has been studied in [13]. Other gen- 
eralizations of the Glauber model consist of, for example, alternating-isotopic 
chains and alternating-bound chains ( [14] for example). 

In [15], an asymmetric generalization of the zero-temperature Glauber model 
on a lattice with boundaries was introduced. There it was shown that in the 
thermodynamic limit, when the lattice becomes infinite, the system shows two 
kinds of phase transitions. One of these is a static phase transition, the other 
a dynamic one. The static phase transition is controlled by the reaction rates, 
and is a discontinuous change of the behavior of the derivative of the station- 
ary magnetization at the end points, with respect to the reaction rates. The 
dynamic phase transition is controlled by the spin flip rates of the particles at 
the end points, and is a discontinuous change of the relaxation time towards the 
stationary configuration. Other phase transitions induced by boundary condi- 
tions have also been studied ( [16-18] for example). Another generalization of 
the Glauber model was introduced in [19]. In this generalization, the processes 
are the same as those of the ordinary Glauber model, but the rates depend on 
three free parameters, rather than one free parameter in the ordinary Glauber 
model. Finally, this model was further generalized to the case where the number 
of interacting sites is more than three and the number of states at each site is 
more than two. This model too violates detailed balance. 

In the present paper an Ising model on a nonuniform lattice with boundaries 
is investigated. Detailed balance is used to propose reaction rates for the system. 
Based on this, the evolution of the expectation values of spins is obtained. The 
time-independent solution to this equation is studied. This solution satisfies a 
homogeneous difference equation of the second order in the bulk, the solution 
to which can be expressed in terms of a transfer matrix. The reactions at the 
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boundaries impose nonhomogeneous (but at most linear) boundary conditions 
on this solution, which could be used to fix the constants appeared in the static 
solution. While it is true that the ordinary Ising model does not exhibit any 
phase transition in finite temperature (the expectation values of the spins van- 
ish if there is no external magnetic field), this is not necessarily the case for 
the model studied here. The expectation values do not vanish as a result of 
inhomogencous boundary conditions. It is shown that in the thermodynamic 
limit (when the size of the lattice tends to infinity) different phases could occur 
for this system, according to whether the eigenvalues of the transfer matrix are 
less than or larger than one. While this transition (or even the existence of 
nonzero static solutions) is invoked by the presence of inhomogeneous bound- 
ary conditions, the detailed form of the boundary conditions affects only the 
coefficients of the eigenvectors of the transition matrix in the static solution. 
So the detailed form of the boundary conditions do not affect the static phase 
portrait of the system. A closed form is obtained for this transfer matrix, and 
some examples are discussed, specially one example in which a phase transition 
is seen. 

The scheme of the paper is as follows. In section 2, the model is intro- 
duced, the rates are determined using the detailed balance criterion, and the 
evolution equation for the spin expectation values is obtained. In section 3, the 
time-independent solution is studied, and the corresponding phase portrait is 
investigated. In section 4, some examples are studied in more detail, specially 
one example which exhibits a static phase transition. Section 5 is devoted to 
the concluding remarks. 

2 One-dimensional Ising model with nonuniform 
coupling constants 

Consider a one-dimensional lattice with (L + 1) sites, labeled from to L. At 
each site, there is a spin variable, Si, which could be +1 for spin up (f), or 
— 1 for spin down (j). These spins in the bulk (sj's with < i < L) interact 
according to the Ising Hamiltonian, 

7~i = ^ ^ Ja 5 a-/i Sa+iii (1) 
a 

where J a is the coupling constant in the link a, and 

,:=\- (2) 

The link a links the sites a — [i and a. + n, so that a ± [i are integers, and a 
runs from fi up to (L — fi) . 

The usual Glauber model gives the dynamics of the Ising model with uniform 
coupling constants ( J a independent of a, and denoted by J) such that the rate 
of a spin flip is determined through its interaction with its two neighboring sites 
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and a heat bath at temperature T. A spin is flipped with the following rates. 



T T T -> T i T and li I — 1 1 i with rate 1 - tanh(2 K), 
T I t -» T T T and | t I — U I with rate 1 + tanh(2 if), 
TU^m and m^ITT with rate 1, (3) 

where 

K := (4) 

and k-Q is the Boltzmann's constant. As it is seen, similar to the Ising model, 
the Glauber model has left-right and up-down symmetries. The Glauber model 
has also a particle reaction-diffusion interpretation. One considers a link with 
different spins at its sites (a domain wall) a particle (•), and a link with same 
spins at its sites (no domain wall) a hole (o). Then the Glauber model turns 
into a reaction-diffusion model: 

o o — » • • with rate 1 — tanh(2 K), 

• • — ¥ o o with rate 1 + tanh(2 K), 

• o ^ o • with rate 1. (5) 

Consider the general case where the coupling constant is not uniform (and the 
interaction is not necessarily nearest neighbor) . Assuming that in each step only 
one spin flips, detailed balance gives 



uj{- ■ ■ , Sj , ► • • • , -sj, ■■■) _ ex p(- • ■ + Ei^j K a s i (~ s j) + ■■■) 

Mil 



£<;(••• >---,Sj,---) cxp(- •• + i& K ij Si Sj + ■■■ ) 

exp(— hj Sj) 



exp(hj Sj) 

where u> is the rate, Kij is defined like (fj| but with Jij (the coupling between 
sites i and j) instead of J, and 

hj:=Y,KijSi- (7) 

As the value of Si is either 1 or —1, any function of Si is at most linear in S{. 
One then arrives at 

exp(/i 3 Sj) = coshhj + Sj sinh/ij. (8) 

Using these, (O gives 

ui(- ■ ■ , Sj, • , — Sj, ■ ■ ■) 1 — Sj tanh hj 



uj{- ■ ■ , — Sj Sj, ■ ■ ■ ) 1 + Sj tanh h 



(9) 



or 



u>{- ■ ■ , Sj, ••■—>•■ ■ , — Sj, ■ ■ •) = Tj (1 — Sj tanh/ij), (10) 
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where Tj-'s are constants. In the simple case of nearest neighbor interaction, one 
has 

Jij = Ji+n — l + Ji-fj, SiJ + l, (11) 

so that ITU1) becomes 



w(- ■ ■ , Sj, ► • • • , -Sj, ■ ■ ■ ) = Tj [1 - Sj tanh(if j _ A( Sj -_i + K j+fl s j+1 )], 

,(12) 

So the spin at the site j flips according to following interactions with the indi- 
cated rates, 

TTT^TIT and Hl^lt-L with rate 1 - tanh^^ + K j+I1 ), 
T I T -> T T T and 1 1 4 -»■ 4 1 -L with rate 1 + tanh(X i _ A1 + K j+(t ), 
T T 4 -» T 4 4 and | 4 T -» 4 T T with rate 1 - tanh(X J -_ AI - K j+It ), 
4TT^44T and T44^TT4 with rate l + tanli(X J -_ M -JC,- + p), 

(13) 

where T^'s have been taken independent of j, and set to one by rescaling the 
time. As it could be expected, the left-right symmetry is violated, but the 
up-down symmetry is not. The particle reaction-diffusion picture turns into fol- 
lowing reaction-diffusion model, which is not left-right symmetric, as expected. 

o o — > • • with rate 1 — tanh(ifj_ Jlt + Kj+p), 

• • — > o o with rate 1 + taiah(Kj-^ + Kj+n), 
o • — > • o with rate 1 — tanJ^Jfj—^ — Kj + ^) 1 

• o — *• o • with rate 1 + tanh(i^j_ p — Kj + ^). (14) 

The evolution equation for the expectation values of the spins in the bulk 
turns out to be 

(sj) = - 2 (sj) + [tarihiKj-n + K j+Ii ) + tanh(i^_ AI - K j+Ii )\ (s^i) 

+ [tanh^-p + K j+lx ) - tanh(^_ M - K j+Il )] (s j+1 ), < j < L. 

(15) 

One has to add two other equations governing the evolution of sq and sl ■ These 
are of the form 

(s ) = a_i + a {.s } + ai 

(s L ) = a L+1 + a L (s L ) + a L _i (.s L -i), (16) 
where a^'s are constants. 
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3 The static solution 

For the static solution ((s) st ), the left hand side of (| 1 5|) vanishes and one obtains 

sinh(2 Kj-fj,) 

\Sj+i)st - — . rjT, ts \ \ s j-i/st 

sinh(2 Kj+n) 

^{ Si ) st , 0<j<Z, (17) 



sinh(2i\j +Al ) 
which can be written as following matrix form 



X j+I1 — Dj Xj-n, 



where 



and 



X a :— 



1 
sinh(2 Kj-^) cosh(2 Kj_^) + cosh(2 Kj +fJ ) 



V smh{2K j+h ) 



sinh(2 K 



One can write Dj as 
where 

and 

Using (HHJ), one arrives at 
where 



I cosh K a sinh i-T Q N 
ysinhXc cosh K a J ' 



sinh K 



3— jj 



cosh ifj. 



\ 



cosh if. 
sinh Kj +fl / 



X a — -Dq,/3 -X"/3, 



A 



so that, 



/ sinh Jf , 
cosh i^Q, 



/3 



M 



cosh if/3 
sinhi^/ 



(18) 
(19) 

(20) 

(21) 
(22) 

(23) 

(24) 

(25) 
(26) 

(27) 
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and 



D a ,p — 



where 



^M Bl/ j - TT-] sinhKp ooshify Kp - M a . p sinh 2 

Aq.^ - -r~— ) sinhifg coshifa co& ^ n __g _ Aq ^ ginll 2 ^ 



(28) 



A a ./3 := tanhii" Q • ■ • taahKp+i, (29) 
M a ^ := tanhi^ Q _i • • ■ tanhii^+i. (30) 



The boundary conditions (|T6|) are 

^■/j = —CL-i, 

Al-^Xl-h = — az+i, (31) 

where 

A a := (a Q _ M a a+(tl ) . (32) 

The steady state profile near the end-site is determined by the eigenvalues 
of the matrix D a fl , where a is some site far from the ends. One has 

v ya r 

where f a is the eigenvector of D a ^ corresponding to the eigenvalue A°, and 
X°'s and XjJ's are the coefficients of expansions of X a and in terms of the 
eigenvectors. It is seen that 

Xl = \ a X«. (34) 

While the exact values of these coefficients are determined by the boundary con- 
ditions pip , one can say whether in the thermodynamic limit these coefficients 
vanish or not, without referring to the exact from of the boundary conditions. 

In the thermodynamic limit, corresponding to each eigenvalue, two cases 
may occur. 

i The eigenvalue A a tends to infinity. In this case 1° tends to zero. 

ii The eigenvalue A a tends to zero or a finite number. In this case 1° is 
generally nonzero. 

Obviously, similar cases occur at the other end point. It is seen that this be- 
havior (some eigenvectors missing or not in the solution corresponding to the 
end points) at one of the end points is independent of the analog behavior at 
the other end. 
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4 Some examples 



Consider some special cases. 

1 Constant coupling: 

K a = K. (35) 

In this case the relation of D a ^ with A ai| g is a similarity transformation. 
So the eigenvalues D at are the diagonal elements of A ai| g, which are 

A 1 = tanh""" K, 

A 2 = coth"^ K, (36) 

showing that one of the eigenvalues is larger and the other is smaller than 
one. So only one of the eigenvectors enters X^. This is regardless of the 
value of K. So there is no static phase transition, as it was seen in the 
case of ordinary (symmetric) Glauber model, [15]. 

2 Periodic coupling: 

K a + m = K a . (37) 

In this case the behavior of the eigenvalues of D a ^ is determined by the 
eigenvalues of D a+m ^ a , which are K a+m , a and A~ +m a . This shows that 
one of the eigenvalues is greater than one and the other is less than one. 
So the situation is similar to that of constant coupling. 

3 Defects in the lattice: No new phenomena is seen, as long as the defects 
are localized, i.e. they are for from the end points. So if there is a lattice 
that has some defects but otherwise is uniform, the static behavior near 
the end points is similar to that of a uniform lattice. 

4 A lattice with different behavior at different end points: The behaviors 
of the static solution near the two ends are independent of each other, 
provided the behavior change occurs far from the end points. So all the 
phenomena seen in previous special cases can be seen at each end point, 
independent of the other end point. 

5 A lattice with different signs of coupling constants: One can define the 
new variables s' and couplings K as 

4 := [sgn(^)...sgn(i^)] Sj , (38) 

K := \K a \ (39) 

where sgn is the sign function, and investigate the system in terms of 
these. So nothing new happens. 

6 A lattice with increasing coupling constants (increasing towards an end 
point): Suppose that K a is an increasing function of a, so that 

lim K a = oo, (40) 
lim Aaj = A, (41) 

a— »oo 
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where A is neither zero nor infinity. Using (|29p , it is seen that (|4"Tj) implies 



lim tanh K a = 1, 

a — >oo 



(42) 



which is equivalent to (|4"D|) . The criterion (|4"Tj) itself, is equivalent to 

lim [cxp(-2 Kr+x) + ■■■+ cxp(-2 K a )] = I, (43) 

a— >oo 

where £ is finite. Assuming that this is the case, it is seen that tends 
to a finite matrix D, where 



D 



A 



A- 



sinh Kp cosh Kg — — - — - — A sinh 2 Kp^ 



sinh Kp cosh Kp — - — A sinh 2 Kp 



A 



(44) 



So in this case both of the eigenvectors remain in X„. 

The last case shows that there is a static phase transition. In one phase only one 
eigenvector remains in X^, while in the other phase both eigenvectors remain. 



5 Concluding remarks 

An Ising model with nonuniform coupling constants on a one-dimensional lat- 
tice with boundaries was studied. Based on detailed balance, the evolution of 
this model was investigated, from which an equation for the static solution was 
obtained. While it is known that ordinary one-dimensional Ising model does not 
exhibit any phase transition at nonzero temperatures, here the inhomogeneity 
at the boundary conditions does permit such a phase transition, while the exact 
form of the boundary conditions do not enter the static phase transition studied 
here. The static phase picture of the system was studied in general, and some 
examples were given including one exhibiting a static phase transition. 

Acknowledgement: This work was partially supported by the research council 
of the Alzahra University. 



8 



References 

[1] R. J. Glauber; J. Math. Phys. 4 (1963) 294. 

[2] K. Kawasaki; Phys. Rev. 145 (1966) 224. 

[3] J. G. Amar & F. Family; Phys. Rev. A41 (1990) 3258. 

[4] T. Vojta; Phys. Rev E55 (1997) 5157. 

[5] R. B. Stinchcombc, J. E. Santos, & M. D. Grynberg; J. Phys. A31 (1998) 
541. 

[6] C. Godrcchc & J. M. Luck; J. Phys. A33 (2000) 1151. 

[7] M. Droz, Z. Racz, & J. Schmidt; Phys. Rev. A39 (1989) 2141. 

[8] B. C. S. Grandi & W. Figucircdo; Phys. Rev. E53 (1996) 5484. 

[9] S. Artz & S. Trimper; Int. J. Mod. Phys. B12 (1998) 2385. 
[10] B. Schmittmann & F. Schmiiscr; Phys. Rev. E66 (2002) 046130. 
[11] B. Schmittmann & F. Schmiiser; J. Phys. A35 (2002) 2569. 
[12] M. Mobilia, R. K. P. Zia, & B. Schmittmann; J. Phys. A37 (2004) L407. 
[13] C. Chatelain; J. Phys. A36 (2003) 10739. 

[14] M. Droz, J. Kamphorst Leal da Silva, A. Malaspinas, & A. L. Stella; J. 
Phys. A20 (1987) L387. 

[15] M. Khorrami & A. Aghamohammadi; Phys. Rev. E63 (2001) 042102. 

[16] M. Henkel & G. Schiitz; Physica A206 (1994) 187. 

[17] M. J. E. Richardson & Y. Kafri; Phys. Rev. E59 (1999) R4725. 

[18] A. Aghamohammadi & M. Khorrami; J. Phys. A34 (2001) 7431. 

[19] A. Shariati, A. Aghamohammadi, & M. Khorrami; Phys. Rev. E64 (2001) 
066102. 



9 



